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Abstract 
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The retarded long-range interaction between neutral systems was first considered by 
Casimir and Polder in their pioneering work in Q|. At large distances atoms interact pre- 



There were various extensions of their result. First of all, Feinberg and Sucher in j3] (see also 
a longer review in P]) expressed the two-photon exchange interaction (dispersion forces) in 
terms of the well-defined physical quantities, the invariant amplitudes which describe the 
elastic scattering of a photon by an atom. These amplitudes are often called dynamic po- 
larizabilities. In this way the authors obtained a formally exact result, as all relativistic 
corrections are included in the dynamic polarizabilities. In the nonrelativistic limit their 
result reduces to that of Casimir and Polder. In the limit of large distances their result is 
expressible in terms of a static electric and magnetic polarizabilities only, similarly to the 
nonrelativistic Casimir-Polder interaction. The calculation of relativistic corrections to the 
dynamic polarizability is not a simple problem. Therefore in later works a different approach 
to dispersion forces was developed which based on a reformulation of nonrelativistic quan- 
tum electrodynamics. Apart from the rederivation of the Casimir-Polder result, there were 
derived closed formulae for higher multipole interactions such as electric quadrupole jj, 0, 0] 
and electric octupole |l7|] for isotropic systems, magnetic dipole and diamagnetic couplings for 
chiral molecules [5|. In this work, using a different reformulation of nonrelativistic quantum 
electrodynamics [8|, we present a systematic derivation of all O(a^) contributions to the 
Casimir-Polder potential of closed shell atoms, including corrections which have not been 
considered so far. The obtained result is expressible in terms of various corrections to the 
nonrelativistic dynamic polarizability, in agreement with the general result of Feinberg and 
Sucher in P]. Moreover, we present the relation of relativistic corrections to the Casimir- 
Polder interaction to the short range nonrelativistic expansion and construct an interaction 
potential which is valid in the whole range of inter-atomic distances. 

Let us first define the a expansion of the interaction potential in the clamped nuclei 
approximation. We use natural units, where c = /i = eo = 1 and denote by m the electron 
mass. The total energy of a system consisting of two neutral atoms is a function of the fine 
structure constant a and the distance R between these atoms 



The nonrelativistic expansion in a depends on the magnitude of R. According to quantum 
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electrodynamics, if R is of the order of an atomic size i? ~ l/(ma), then this expansion at 
constant ma R takes the form 



E{a, mR) = E^^) (maR) + E^^^ {maR)+ E^^^ ( ma R) + 0{ 
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where E^"^^ is the nonrelativistic energy of order m of a systems of two atoms including 
Coulomb interactions between all electrons, E^^^ is the leading relativistic correction of order 
m a'^, which is given by the Breit-Pauli Hamiltonian 6H P|. We include below only the terms 
which do not vanish for closed shell atoms 
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where the sum goes over all electrons of both atoms. -E*^^^ is the QED correction of order 



ma^. It consists of various terms, among others the Araki-Sucher term 
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is dominating at large atomic distances, 
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On the other hand, if R is of the order of atomic transition wavelength, namely R ~ 
l/{ma^) then the a expansion at constant ma^ R takes a completely different form 

E{a, mR) = Efj.f,f,{a) + Ecp{mL a ^R) + 5^^^Ecpima^R) + ... (6) 

where Efj-ee is the energy of separate atoms, Ecp is a Casimir- Polder potential jl| and S^'^^'Ecp 
is the leading relativistic correction, which is the subject of this work. The form of this 
expansion results from the long wavelength formulation of quantum electrodynamics, see 
Ref. fl. This relativistic correction to the interaction energy at large atomic distances is 
obtained from the effective interaction Hj of an atom with the slowly varying electromagnetic 
field |S| 



Hi = y -efa - E - - iriri E\: r^rl E\ 
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where L = r x p, E = E{0), B = B{0) are fields at the position of nucleus, and spin 
dependent terms have been neglected as we consider only closed shell atoms. The sum in 



Eq. (|7j) goes over all electrons of one atom. For simplicity, we will assume this sum is 
present implicitly in all the formulas below. The leading Casimir-Polder interaction comes 
from the two-photon exchange with the interaction —efa ■ E. Using the temporal gauge for 
the photon propagator = 0, it is 3] 

EcP = -"^ r ^^<aC^'eb9''9'\ (8) 
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and tu-integration is assumed along the Feynman contour. This integration contour is de- 
formed to imaginary axis by the replacement uj = i\. The fc-integral leads to 
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This result has been obtained by Casimir and Polder in We consider here the 
correction and represent it as a sum of five terms 

5^^^EcP = df^EcP + df^EcP + df^EcP + df^EcP + df^EcP . (13) 

6^^Ecp is due to the Breit-Pauli correction to H, E and state (f) in Eq. (fT^ . For simplicity 
we consider corrections only to the atom A, therefore only one matrix element in Eq. ()12p 

(2) 

is to be modified in 6q Eqp, according to 
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. (i7-E)2 + A2 / \ (E-H)' {H-Ef + X^ 

This con.ec«o. has .ecen.y bee. ec.de.ed in Ref. Q. All ..nai.ing eo^eCo,. .„ 

the Casimir-Polder interaction energy are obtained by modification of a dipole interaction 
—ef-E by various couplings as given by Eq. 0. (5f ^EcP comes from the quadrupole term 
-e/2 (rl ri - 5'^ r^/S) E^ in Eq. Q, 

^ """^ 120 io 27r\ {H-Ef + X^ I a\ {H - Ef + Ib 

where 
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After contracting i,j, k indices it becomes 

6^'^Er. - r dxlr-- 2{H-E) \ I 2{H - E) 

SOvrio ^^V {H-Ef + X^' I a\ {H - Ef + X^' I b 
^ -2Xr[^ J_ 27 84 162 180 90 \ 

This result was first obtained by Jenkins, Salam and Thirunamachandran in [5]. Ecp is 
the correction that comes from the term — e/30 r* E'^^^ in Eq. (|7j). It is very similar to the 
dipole-dipole interaction and can easily be obtained on the basis of Eq. (fT^ 

^&^-2\R , 2 5 6 3 \ . 
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^3 Ecp is due to — e/ (6 m) {L\ + L\) and can be regarded as another correction to 
the electric dipole coupling 

,pEc, = -1^ r^A/ f,,, „ f-l)(Lx F-FxL)\ ( f"-f, r\ 
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The last correction Sf^Ecp is due to e^/(8m^) (r^ x B)^ and reads 
„2 .00 / H-E \ 
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The complete O(a^) correction is a sum of Eqs. ()14I18I19I20I21|) as given by Eq. fjl3p . Here, 
5i Ecp comes from interaction between the electric dipole and the electric quadrupole 

polarizabilities, 6i^^Ecp is the interaction energy of the electric dipole polarizability with 

(2) (2) (2) 

the magnetic susceptibility, and Ecp + Sq Ecp + 62 Ecp + 6l^ Eqp is the interaction energy 
between electric dipole polarizabilities with the relativistic correction 
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We have not found in the literature the complete formula for the leading relativistic correc- 
tion to the electric dipole polarizability of the closed shell atoms as that in the Eq. (j^Hj) . 

Let us now consider the large and small R limit of the interaction energy. At large R, 
the Sq'^Ecp and S^^Ecp contribute to the coefficient, but it is only a small correction 
on the top of Ecp- Much more interesting is a small R expansion of Ecp and 5^'^^Ecp 
and its relation to the large R expansion of energy as a function of a and ma R. This 
relation has been first considered by Meath and Hirschfelder in 3|- The large R expansion 
of E^^\ma R) from Eq. Q reads 
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where Cp'* are dimensionless constants. The relation to the small R expansion of 
Ecp{m a ^R) and S^^^Ecp{ma^ R) from Eq. (jH} is the following 
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Since both expansions of E(a,mR) involve the same coefficients C-'\ one can write the 



general formula 



E{a,mR) = £^free(a) — ^ma* 
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We have checked this by equivalence of Cp^ coefficients as obtained from these two different 

expansions, and they are equal to (in atomic units) 
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where Na and Nb are the number of electrons in the atom A and B respectively. If S^'^^Ecp 
includes contributions from the atom B, then coefficients Cg^\ Cg^^ and Cg^^ should include 
corresponding terms obtained by the replacement A S. 

The Cj'^ coefficients allow one to obtain a convenient form of the interaction potential 
in the whole region of the atomic distance R, as long as these atoms do not overlap. The 
minimal version of this potential is 



maR) + Ecp{ma R) + ma 
and the most accurate version using present result is 
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In summary, the purpose of this work was the derivation of a complete cP' correction to 
the Casimir-Polder potential in order to obtain a more accurate description of inter-atomic 
interactions in the region where the electron wave functions from different atoms do not 
overlap. The obtained result can be used for the precise calculation of the scattering length 
and highly excited vibrational levels of light molecules. Particularly interesting is the helium 



ny II 



dimer which existence has been confirmed as recently as in 1994 [16[. Its dissociation energy 
has an extremely small value of 1 mK, while the mean internuclear distance is as large as 50 
A. Its existence can be associated to the long range attraction between the helium monomers. 
Since the minor perturbations of the interaction potential result in significant changes in the 
description of the nuclear motion, the potential in the large range of inter-atomic distances 
with an accuracy of the order of 1 mK is needed, which is the magnitude of relativistic and 
QED effects [17]. 
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